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Association of self-reported race with AIDS death in 
continuous HAART users in a cohort of HIV-infected 

women in the United States 

Kerry Murphy 3 , Donald R. Hoover b , Qiuhu Shi c , Mardge Cohen d,e , 
Monica Gandhi f , Elizabeth T. Golub 8 , Deborah R. Gustafson h '', 
Celeste Leigh Pearce', Mary Young k and Kathryn Anastos a 

Objective: To assess the association of race with clinical outcomes in HIV-positive 
women on continuous HAART. 

Design: Prospective study that enrolled women from 1 994 to 1 995 and 2001 to 2002. 

Setting: Women's Interagency HIV Study, a community-based cohort in five US cities. 

Participants: One thousand, four hundred and seventy-one HIV-positive continuous 
HAART users. 

Main outcome measures: Times to AIDS and non-AIDS death and incident AIDS- 
defining illness (ADI) after HAART initiation. 

Results: In adjusted analyses, black vs. white women had higher rates of AIDS death 
[adjusted hazard ratio (aHR) 2.14, 95% confidence interval (CI) 1.30, 3.50; P = 0.003] 
and incident ADI (aHR 1.58, 95% CI 1.08, 2.32; P=0.02), but not non-AIDS death 
(aHR 0.91 , 95% CI 0.59, 1 .39; P = 0.65). Cumulative AIDS death incidence at 1 0 years 
was 17.3 and 8.3% for black and white women, respectively. Other significant 
independent pre-HAART predictors of AIDS death included peak viral load (aHR 
1.70 per log 10 , 95% CI 1.34, 2.16; P<0.001), nadir CD4+ cell count (aHR 0.65 per 
100 cells/|xl, 95% CI 0.56, 0.76; P< 0.001), depressive symptoms by Center for 
Epidemiology Studies Depression score at least 16 (aHR 2.10, 95% CI 1.51, 2.92; 
P < 0.001 ), hepatitis C virus infection (aHR 1 .57, 95% CI 1 .02, 2.40; P= 0.04), and HIV 
acquisition via transfusion (aHR 2.33, 95% CI 1.21, 4.49; P=0.01). In models with 
time-updated HAART adherence, association of race with AIDS death remained 
statistically significant (aHR 3.09, 95% CI 1.38, 6.93; P = 0.006). 

Conclusion: In continuous HAART-using women, black women more rapidly died 
from AIDS or experienced incident ADI than their white counterparts after adjusting for 
confounders. Future studies examining behavioral and biologic factors in these women 
may further the understanding of HAART prognosis. 
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Introduction 



Globally, women represent more than 50% of persons 
living with HIV infection (PLWH), and HIV/ AIDS is the 
leading cause of death worldwide for women of 
reproductive age [1,2]. In the United States, African 
American women are 19-fold more likely than white 
women to be diagnosed with HIV [3], and 2007 United 
States Centers for Disease Control and Prevention data 
ranked HIV in the top 10 causes of death for African 
American and Hispanic women aged 20—54 years [4,5]. 
Declines in AIDS-related mortality have been less steep in 
black women than in most other groups [6,7]. Many 
studies of safety and efficacy of antiretroviral therapy 
(ART) have been conducted in men of European descent; 
it is unclear if there is a survival difference by race in 
PLWH treated with HAART. Whereas some studies have 
found no differences by ancestry in risk of death or AIDS- 
defming illness (ADI) in HAART users [8-10], others 
have observed higher rates of all-cause mortality or 
decreased rates of virologic suppression in people of 
African descent, but did not assess or find any difference 
in AIDS-specific mortality or adherence [11 — 16]. 

Postulated explanations for higher AIDS mortality rates 
among blacks, if present, include decreased access to care, 
lower HAART adherence, and socioeconomic, beha- 
vioral, genetic, or other biologic factors [17—20]. In a 
previous study from the Women's Interagency HIV Study 
(WIHS), we found a clinically but not statistically 
significant survival advantage in white compared to 
black women on continuous HAART: adjusted hazard 
ratio (aHR) 0.49 (P=0.19) for death from AIDS [21]. 
We extended the WIHS study into the era of modern 
HAART using follow-up and events from 1995 to 2009, 
hypothesizing that black women would have higher rates 
of AIDS death and ADI compared to white women. 



Methods 
Study population 

Women's Interagency HIV Study is a multicenter 
prospective observational cohort study of 2843 HIV- 
infected and 975 HIV-uninfected women enrolled in 
1994-1995 and 2001-2002. Informed consent was 
obtained from participants, and the Institutional Review 
Boards of the six collaborating institutions approved the 
study. Study participants were interviewed and examined 
semiannually; numerous laboratory tests and studies were 
performed [22]. This analysis included all WIHS 
participants who initiated HAART after WIHS enroll- 
ment and reported HAART use at 70% or more of 
subsequent visits. Although individual women were 
allowed to have up to four consecutive missed visits or 
visits off HAART (i.e. in cases of otherwise long-term 
consistent use), less than 5% of the total semiannual 



person-visits were off HAART. HAART was defined as 
any three-drug antiretroviral regimen, one of which must 
have been a protease inhibitor, a non-nucleoside reverse 
transcriptase inhibitor, one of the nucleoside reverse 
transcriptase inhibitors abacavir or tenofovir, an integrase 
inhibitor (raltegravir), or an entry inhibitor (maraviroc or 
enfuvirtide) [23]. 

Exposure variables 

The primary exposure was self-categorized race: white 
(Hispanic and non-Hispanic), black (Hispanic and non- 
Hispanic), and other, who were predominantly women 
who self-identified as Hispanic, but not white or black. 
Women self-categorized their race and ethnicity by 
responding to interviewer-administered questions: one 
querying whether they were Hispanic or not, and another 
asking if they were black/ African American, white/ 
Caucasian, American Indian/ Alaskan Native, Asian, 
Native Hawaiian/ Other Pacific Islander, or other. Other 
variables examined included date of HAART initiation 
defined as the mid-date between the visit at which 
consecutive HAART use was first reported and the 
previous visit; age at HAART initiation; time of 
enrollment in WIHS (1994-1995 or 2001-2002); risk 
category for HIV acquisition [injection drug use (IDU), 
heterosexual exposure, transfusion]; and pre-HAART 
initiation variables of illicit drug use, smoking tobacco, 
both in the past 6 months, annual income $12 000 or less, 
having graduated from high school, self-reported prior 
ADI; hepatitis C virus (HCV) infection defined as 
positive antibody with viremia; nadir CD4 + cell count 
and peak log 10 HIV-1 RNA, which reflected the lowest 
CD4 + cell count; and highest HIV viral load measured at 
any pre-HAART visit. We also included pre-HAART 
clinical depressive symptom burden defined as Center for 
Epidemiology Studies Depression score (CES-D) at least 
16. As others have done [24—27], we divided the 20 CES- 
D items into depression component subscores: somatic, 
positive, negative, and interpersonal. Each of these 
subscores was dichotomized about the closest value to 
the sample median to obtain somatic, negative, positive, 
and interpersonal depression subscales of at least 7, 5, 4, 
and 1, respectively. We adjusted for WIHS site in all 
analyses. Starting from October 2002, WIHS included 
measures for assessment of self-reported adherence 
(defined as >95% use of all prescribed antiretrovirals) , 
which were included in a separate analysis. 

Outcome variables 

Primary outcomes were time from date of first visit after 
HAART initiation to AIDS death, non-AIDS death, and 
incident ADI. Ascertainment and classification of deaths 
in WIHS have been described previously [28,29]. For this 
analysis, deaths from infection were included as AIDS 
deaths. ADIs were self-reported and classified as incident 
as previously described [30]. To account for loss to 
follow-up, women without study outcomes were 
censored at the earlier of 2 years after their last HAART 
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use visit or 31 December 2009 as ascertainment of death 
beyond that date may not have been complete. We 
analyzed self-reported HAART adherence at each 6- 
month visit. 

Laboratory methods 

Plasma HIV-1 RNA was measured by isothermal nucleic 
acid sequence-based amplification (NASBA/Nuclisens; 
Organon Teknika Corp., Durham, North Carolina, 
USA) with a lower limit of detection of 80 copies/ ml. 
Lymphocyte subsets were quantified using standard flow 
cytometric methods in laboratories participating in the 
NIH/NIAID Flow Cytometry Quality Assessment 
Program [31]. HCV RNA was measured on frozen 
specimens from HCV antibody-positive women using 
either the COBAS Amplicor Monitor 2.0 assay or the 
COBAS Taqman assay (Roche Diagnostics, Branchburg, 
New Jersey, USA) [32J. 

Statistical analysis 

Racial groups were compared by means or proportions of 
potential confounding covariates. We performed unad- 
justed competing risk Kaplan— Meier and adjusted 
proportional hazards survival models of time to AIDS 
death, non-AIDS death, and incident ADI by racial group 
[33]. In models restricted to 1255 women with adherence 
data, adherence to HAART (available from October 2002 
forward) was included as a time-dependent variable. 
Women without study outcomes were censored at the 
earlier of 2 years after their last HAART use visit or 31 
December 2009. Multivariate proportional hazards 
analyses were performed with backwards stepwise 
selection among the following variables assessed at the 
pre-HAART visit: race, age, HIV exposure category, 
income, education, date of WIHS enrollment, WIHS 
site, year of HAART initiation, illicit drug use, smoking 
status, depressive symptoms, nadir CD4 + , peak HIV-1 
RNA, pre-HAART ADI, HCV status, and adherence. As 
the primary exposure of interest, race was forced into the 
stepwise models, as were age and smoking tobacco. 
Logistic regression models of adherence from October 
2002 onwards were fit using the same predictor variables, 
and generalized estimation equations (GEEs) to adjust for 
within-person collinearity of repeated visits [34] . Multi- 
variate models of adherence were fit using the same 
stepwise approach (without adherence as a predictor) 
described above. A sensitivity analysis assessing collinear- 
ity of HCV and IDU was also performed. 



Results 



Of the 2090 HIV-positive WIHS participants who 
reported HAART use, 1471 met inclusion criteria for 
continuous HAART use: 823 (55.9%) were black, 334 
(22.7%) white, and 314 (21.3%) reported their race as 
other. During the study follow-up, 167 women died of 



AIDS, 136 died of non-AIDS causes, and 11 with 
unknown cause of death. Most deaths of unknown cause 
were recent without enough information available to 
allow adjudication for cause. Clinical and sociodemo- 
graphic characteristics of the study population are shown 
in Table 1 . Black compared to white women were older at 
HAART initiation, more likely to have used illicit drugs, 
have HCV, and smoke tobacco. Women reporting their 
race as other initiated HAART earlier in calendar time 
and had lower income and educational attainment. There 
was no statistical difference in the prevalence of depressive 
symptoms at the pre-HAART visit across racial groups. 
HAART adherence was significantly lower in 
black women. 

Cumulative incidence of death 

Figures 1 and 2 show cumulative incidence of AIDS and 
non-AIDS death, respectively, stratified by race. Ten years 
after HAART initiation, black women had the highest 
cumulative incidence of AIDS death, 17.3 vs. 8.3% for 
whites and 14.2% for others (Fig. 1), whereas cumulative 
incidence of non-AIDS death was similar in blacks and 
whites (13.1 and 12.3%, respectively) and lower (7.7%) in 
women reporting their race as other (Fig. 2). 

Associations with death from AIDS 

In univariate and multivariate analyses, black women 
experienced more rapid death from AIDS than white 
women. In multivariate models not adjusting for 
adherence (middle column, Table 2) the aHR for AIDS 
death in black vs. white women was 2.14 (95% CI 1.30, 
3.50; P= 0.003). Other significant independent pre- 
dictors of AIDS death were higher peak pre-HAART 
viral load (aHR 1.70 per log 10 HIV RNA increase, 95% 
CI 1.34, 2.16; P< 0.001), HCV status (aHR 1.57, 95% 
CI 1.02, 2.40; P=0.04), and depressive symptoms (aHR 
2.10, 95% CI 1.51, 2.92; P< 0.001). As shown in 
univariate models (left column, Table 2), negative, 
positive, interpersonal, and somatic depressive symptoms 
qualitatively had equivalent associations with time to 
AIDS death; thus only CES-D at least 16 was included in 
the multivariate model because of issues of collinearity 
and multiple comparisons. Women who acquired HIV 
via transfusion vs. heterosexual exposure were more likely 
to die from AIDS (aHR 2.33, 95% CI 1.21, 4.49; 
P=0.01). Higher pre-HAART nadir CD4 + cell count 
(aHR 0.65 per 100 cells/jxl increment, 95% CI 0.56, 
0.76; P< 0.001) and later year of HAART initiation 
(aHR 0.91 per year, 95% CI 0.85, 0.97; P= 0.005) 
protected against AIDS death. 

In both univariate and multivariate subanalyses of the 
1255 women with adherence data, at least 95% self- 
reported adherence to HAART strongly protected against 
AIDS death (aHR 0.33, 95% CI 0.21, 0.53; P< 0.001) 
(right column, Table 2). When time-updated adherence 
was included in the model, the association of African 
American race with AIDS death persisted (aHR 3.09, 
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Table 1. Participant characteristics in percentages (%) or means (SDs). 
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CES-D, Center for Epidemiology Studies Depression score; HCV, hepatitis C virus; IDU, injection drug use; WIHS, Women's Interagency HIV 
Study. 

a AII person-visits from October 2002 onwards for which adherence was self-reported in 1255 women. P-values are from generalized estimation 
models. 

b See Table 2 for statistical comparison of times to AIDS and non-AIDS deaths. 



Cumulative incidence of AIDS death over time 
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Fig. 1. Cumulative incidence of AIDS death over time. Lines 
correspond to the cumulative incidence of AIDS death stra- 
tified by race. The red line represents black women, the blue 
line represents white women and the green line represents 
women who reported their race as other. 



Cumulative incidence of non-AIDS death over time 
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Fig. 2. Cumulative incidence of non-AIDS death over time. 

Lines correspond to the cumulative incidence of non-AIDS 
death stratified by race. The red line represents black women, 
the blue line represents white women and the green line 
represents women who reported their race as other. 
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95% CI 1.38, 6.93; P= 0.006). Adjustment for adherence 
did not statistically or qualitatively attenuate the pre- 
viously described association of other variables with time 
to AIDS death. 

Associations with non-AIDS death 

There was no significant association of race with non- 
AIDS death; aHR 0.91, 95% CI 0.59, 1.39; P= 0.65 and 
aHR 0.75, 95% CI 0.42, 1.34; P=0.34 for black and 
other vs. white, respectively (right column, Table 3) in 
multivariate analyses. Significant independent predictors 
of the incidence of non-AIDS death were: HCV status 
(aHR 2.45, 95% CI 1.67, 3.60; P< 0.001), older age at 
HAART initiation (aHR 1.19 per 5 years, 95% CI 1.06, 
1.33; P= 0.004), smoking (aHR 1.61, 95% CI 1.06, 2.44; 
P= 0.03), and annual income $12 000 or less (aHR 1.57, 
95% CI 1.05, 2.35; P— 0.03). There was no association of 
HAART adherence or depressive symptoms with non- 
AIDS death in univariate or multivariate models, 
respectively; thus the multivariate model of time to 



non-AIDS death adjusting for time-dependent adherence 
to HAART is not presented. 

Incident AIDS-defining illness 

The proportion of women who reported an incident ADI 
at 10 years was 21.3, 13.9, and 16.6% for black, white, and 
other women, respectively (see Supplemental Digital 
Content 1 , Figure of Cumulative Incidence of Incident 
AIDS-Defming Illness, http://links.lww.com/QAD/ 
A386). In multivariate analyses, black women (aHR 
1.58, 95% CI 1.08, 2.32; P=0.02) and women who 
reported their race as other (aHR 1.60, 95% CI 1.01, 
2.54; P— 0.045) were more likely than white women to 
report an incident ADI. Incident ADIs were also more 
common in women with a higher pre-HAART peak 
viral load (aHR 1.51 per log 10 increment, 95% CI 1.26, 
1.81; P< 0.001), prior history of ADI (aHR 2.19, 95% 
CI 1.63, 2.96; P< 0.001), and current smokers (aHR 
1.65, 95% CI 1.24, 2.20; P< 0.001), and less common in 
women who initiated HAART in more recent 



Table 3. Predictors of non-AIDS death hazards ratios (HRs) or adjusted hazard ratios (aHRs). a 

Univariate model (N=1471) Multivariate model b ' c (N=1471) 



HR (95% CI) P-value aHR (95% CI) P-value 



Demographics 
Race 

Black vs. white 1.25(0.81,1.91) 0.31 0.91(0.59,1.39) 0.65 

Other vs. white 0.79(0.45,1.40) 0.42 0.75(0.42,1.34) 0.34 

Age at visit per 5 years 1.24(1.12,1.36) <0.001 1.19(1.06,1.33) 0.004 

Mode of HIV acquisition 

IDU vs. heterosexual 3.15(2.18,4.55) <0.001 

Transfusion vs. heterosexual 2.66(1.14,6.22) 0.02 

Pre-HAART annual income <$1 2 000 2.20(1.50,3.21) <0.001 1.57(1.05,2.35) 0.03 

Education: high school graduation or beyond 0.61 (0.43, 0.86) 0.005 

WIHS year of enrollment (1 994-1 995 0.48 (0.23, 0.98) 0.044 

vs. 2001-2002) 

Calendar year of HAART initiated 0.95 (0.90, 1 .01 ) 0.1 0 

Behaviors pre-HAART 

Pre-HAART illicit drug use 1 .82 (1 .20, 2.78) 0.005 

Pre-HAART smoking 2.63(1.81,3.84) <0.001 1.61(1.06,2.44) 0.03 

Pre-HAART depressive symptoms 1.73(1.23,2.44) 0.002 1.40(0.98,2.00) 0.06 

CES-D 16+ vs. <16 

Pre-HAART CES-D subscale depression 

Negative subscale <4 vs. >4 1.65(1.17,2.33) 0.004 NA d NA d 

Somatic subscale <6 vs. >6 1.60(1.14,2.25) 0.007 NA d NA d 

Positive subscale <3 vs. >3 1.75(1.24,2.46) 0.001 NA d NA d 

Interpersonal subscale 0 vs. >0 1.46(1.04,2.05) 0.03 NA d NA d 

Clinical characteristics pre-HAART 

Pre-HAART nadir CD4+ per 100 cells/(J 0.88 (0.78, 1.00) 0.051 

Pre-HAART log, 0 peak VL 1.17 (0.94, 1 .45) 0.1 6 

Pre-HAART AIDS-defining illness 1.68(1.18,2.38) 0.004 

Pre-HAART HCV infection 3.58(2.51,5.11) <0.001 2.45(1.67,3.60) <0.001 

Post-HAART events/outcomes 

Adherence to HAART >95% vs. <95% b ' e 0.80 (0.49, 1 .29) 0.36 



CES-D, Center for Epidemiology Studies Depression score; CI, confidence interval; HCV, hepatitis C virus; IDU, injection drug use; WIHS, 
Women's Interagency HIV Study; VL, viral load. 
a From proportional hazards model. 

b From stepwise selection proportional hazards models with race, WIHS site, age, and smoking forced in 

c Note: As adherence was not associated with non-AIDS death in unadjusted analysis or multivariate analysis, the multivariate model with 
adherence is not reported here as was done in Table 2. 

d NA - only CES-D and not the specific subscales was considered for the multivariate model. 

Restricted to survival from October 2002 through December 2009 for which self-reported adherence was obtained. 
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calendar time (aHR 0.75 per year, 95% CI 0.70, 
P< 0.001). 

Analyses excluding illicit drug users 

To eliminate possible confounding by drug use, we 
performed analyses excluding women who reported any 
history of using crack, cocaine or heroin, or a history of 
IDU at the pre-HAART visit. With drug users (33.1% of 
participants) excluded, the association of black (vs. white) 
race with AIDS death was larger (aHR 3.10, 95% CI 1.47, 
6.53; P= 0.003); black women again more rapidly 
experienced an incident ADI (aHR 1.90, 95% CI 
1.12, 3.24; P=0.02) in this subset. With drug users 
excluded, there was again no association of race with non- 
AIDS death (aHR for African American race 1.32, 95% 
CI 0.65, 2.66; P=0.43) and little change in the 
association of other predictors with AIDS deaths (see 
Supplemental Digital Content 2, Supplemental Table 1, 
http://links.lww.com/QAD/A386 Predictors of AIDS 
Death with Illicit Drug Users Excluded and Supple- 
mental Digital Content 3, http://links.lww.com/QAD/ 
A386, Supplemental Table 2, http://links.lww.com/ 
QAD/A386 Predictors of Non-AIDS Death with Illicit 
Drug Users Excluded). 

Sensitivity analysis of injection drug use and 
hepatitis C virus 

There was some association between HIV acquisition via 
IDU and HCV infection; 73% of IDU vs. 25% of women 
who acquired HIV via heterosexual exposure were 
infected with HCV. A sensitivity analysis using three 
different models of predictors of AIDS death (HCV 
infection alone, IDU alone, and both) indicated that 
HCV infection and IDU were not collinear in their effect 
on AIDS deaths (data not shown). 

Associations of race and other covariates with 
adherence 

We fit GEE models assessing predictors of adherence at all 
visits at which HAART use was reported. In univariate 
and multivariate models, black women were less likely 
than white women to report at least 95% HAART 
adherence [adjusted odds ratio (aOR) 0.66, 95% CI 0.53, 
0.83; P< 0.001] (Table 4). Women with depressive 
symptoms were also less likely to adhere to HAART 
(aOR 0.77, 95% CI 0.65, 0.92; P = 0.004); each of the 
depression subscales (negative, somatic, positive, and 
interpersonal) was qualitatively associated with non- 
adherence with similar strength. Other independent 
predictors of lower adherence were illicit drug use (aOR 
0.75, 95% CI 0.58, 0.98; P=0.03), HCV status 
(aOR 0.77, 95% CI 0.63, 0.95; P= 0.02), and smoking 
(aOR 0.73, 95% CI 0.61, 0.88; P< 0.001). There was 
no association between duration of HAART use and 
adherence. 



In this large study of continuous HAART users with 
varied ancestry, we found that black women were twice as 
likely as white women to experience adverse HIV clinical 
outcomes, specifically death from AIDS and incident 
ADI, even after adjusting for multiple potential 
confounding characteristics including illicit drug use, 
depressive symptoms, and adherence. These findings of 
strong independent associations of African American race 
with AIDS death and incident ADI after controlling for 
important potential confounders raise the possibility that 
both genetic and biologic factors, or either of them, could 
be substantially influencing HAART efficacy. These 
differences seem specific to HIV as there were no 
differences by race in non-AIDS deaths. 

Collection of self-reported HAART adherence data in 
WIHS after 2002 allowed us to assess the contribution of 
adherence to higher death rates in black women. 
Although black women were less likely to adhere to 
HAART (aOR 0.66) than white women, multivariate 
models including time-updated adherence did not 
mitigate this group's faster rate of death from AIDS. 
Further, in models excluding illicit drug users, the 
association of AIDS deaths with ancestry was even larger 
(hazard ratio 3.09 vs. 2.14 for analyses excluding and 
including drug users, respectively). These findings suggest 
that factors other than modifiable behaviors may be 
contributing to racial disparities in AIDS clinical 
outcomes. 

The finding of a strong independent association of 
African American race with AIDS-specific clinical 
outcomes after accounting for numerous potential 
confounding factors raises the question of what other 
potential mediators may be contributing. There is the 
possibility of unmeasured confounders that we did not 
capture through our chosen covariables. For example, we 
analyzed annual income and education level to assess 
socioeconomic status, but perhaps variables such as 
neighborhood, family structure, and employment status 
would have provided additional measures. Similarly, 
unmeasured behavioral factors may have varied by race 
and could influence mortality. Because all participants by 
definition were accessing healthcare to obtain HAART, 
lack of access to care is substantially obviated as a 
confounder; however, measurement of health insurance 
status (public vs. private) and patient-provider relation- 
ships may be informative. 

The study results suggest that factors other than 
demographics, socioeconomic status, and modifiable 
behaviors may be contributing to racial disparities in 
AIDS clinical outcomes. Multiple sources of evidence 
suggest that genetically determined diversity in immu- 
nogenetic and pharmacokinetic parameters is associated 
with clinical findings in HIV-infected persons using 
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Table 4. Predictors of at least 95% adherence to HAART odds ratios (ORs) or adjusted odds ratios (aORs).' 



Univariate model (N=1255) 



Multivariate model b (N= 1255) 





OR (95% CI) 


P-value 


aOR (95% CI) 


P-value 


Demographics 










Race 










Black vs. white 


0.63 (0.51, 0.80) 


<0.001 


0.66 (0.53, 0.83) 


<0.001 


Other vs. white 


0.80 (0.61, 1.05) 


0.10 


0.76 (0.58, 1.00) 


0.051 


Age at visit per 5 years 


1.06 (1.00, 1.11) 


0.04 


1.10 (1.04, 1.17) 


0.001 


Mode of HIV acquisition 










IDU vs. heterosexual 


0.75 (0.62, 0.90) 


0.002 






Transfusion vs. heterosexual 


1.38 (0.83, 2.29) 


0.22 






Pre-HAART annual income <$12 000 


0.85 (0.71, 1 .01) 


0.06 






Education: high school graduation or beyond 


1.00 (0.84, 1.20) 


0.96 


0.86 (0.73, 1 .02) 


0.09 


WIHS year of enrollment (1 994-1 995 vs. 2001 -2002) 


1.04 (0.84, 1.28) 


0.72 






Calendar year HAART initiated 


1.00 (0.98, 1.02) 


0.96 






Behaviors pre-HAART 










Pre-HAART illicit drug use 


0.57 (0.44, 0.73) 


<0.001 


0.75 (0.58, 0.98) 


0.03 


Pre-HAART smoking 


0.63 (0.53, 0.74) 


<0.001 


0.73 (0.61, 0.88) 


<0.001 


Pre-HAART depressive symptoms CES-D 16+ vs. <16 


0.72 (0.61, 0.86) 


<0.001 


0.77 (0.65, 0.92) 


0.004 


Pre-HAART CES-D subscale depression 










Negative subscale <4 vs. >4 


0.73 (0.61, 0.86) 


<0.001 


NA C 


NA C 


Somatic subscale <6 vs. >6 


0.65 (0.55, 0.77) 


<0.001 


NA C 


NA C 


Positive subscale <3 vs. >3 


0.77 (0.65, 0.91) 


0.002 


NA C 


NA C 


Interpersonal subscale 0 vs. >0 


0.82 (0.69, 0.98) 


0.03 


NA C 


NA C 


Clinical characteristics pre-HAART 










Pre-HAART nadir CD4 + per 100 cells/u,l 


1.08 (1.02, 1.15) 


0.01 


1.07 (1.01, 1.14) 


0.02 


Pre-HAART log 10 peak VL 


0.91 (0.82, 1.01) 


0.07 






Pre-HAART AIDS-defining illness 


0.77 (0.65, 0.91) 


0.002 






Pre-HAART HCV infection 


0.72 (0.60, 0.87) 


<0.001 


0.77 (0.63, 0.95) 


0.02 



CES-D, Center for Epidemiology Studies Depression score; CI, confidence interval; HCV, hepatitis C virus; IDU, injection drug use; WIHS, 
Women's Interagency HIV Study; VL, viral load. 

a Based on GEE logistic regression. Study visit was also considered but was not statistically associated with adherence. 
b From backwards selection GEE models with race, WIHS site, WIHS visit, age, and smoking forced in. 
C NA - only CES-D and not the specific subscales was considered for the multivariate model. 



HAART and may be contributing to our findings of a 
two-fold higher risk in black vs. white women of death 
from AIDS, the single most important clinical outcome. 

Host immunogenetic and pharmacogenetic variability are 
both known to influence clinical outcomes in HIV 
disease. Several studies have shown that specific human 
leukocyte antigen complex (HLA) class 1 alleles, which 
are known to vary by ancestry, are associated with HIV 
disease progression [35-38], including elite suppression 
[39]. Further, recent studies have indicated that some 
alleles (e.g. HLA B*5701 and Bw4), which are protective 
for progression in untreated HIV disease, may have 
deleterious effects in patients taking HAART, with a 
lower likelihood of virologic suppression or CD4 
response [40—43]. Chemokine receptor gene poly- 
morphisms have also been associated with differential 
clinical outcomes on HAART based on ancestry (e.g. 
CCR5 2459G allele) [44]. Pharmacogenetic factors that 
can influence drug metabolism [45-48] and transport 
[49-51] are known to vary by ancestry, but their effects 
on response to HAART are not clear [52—63]. 

We found that higher depressive symptom scores were 
associated with AIDS death after adjusting for other 
covariates. The significant overlap of items of the CES-D 



with somatic symptoms and advanced disease (fatigue, 
anorexia, lower energy, and lack of mental acuity) could 
be confounding the association of depression with poor 
clinical outcomes, as has been found in prior studies [64]. 
However, we found that all components of the CES-D 
score including somatic, negative, positive, and inter- 
personal subscales were significantly associated with AIDS 
death at qualitatively similar levels (left column, Table 2). 
Thus the strong association of depressive symptoms with 
AIDS death did not result just from the somatic 
components within the CES-D score. This finding is 
clinically significant as depression is potentially modifiable 
through treatment [65—67]. 

Although the black women were more likely to smoke 
tobacco, and tobacco use was associated with AIDS death 
in univariate analysis, its use was not associated with AIDS 
death in the multivariate models. Tobacco use was, 
however, associated with non-AIDS death as expected. 

Strengths of our study include its large sample size, 
representativeness of the WIHS cohort, and the long 
follow-up. This is one of the largest studies to date to 
assess clinical outcomes by ancestry and adjusts for 
confounding behaviors and characteristics that are not 
typically measured in clinical settings. The WIHS dataset 
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includes measurements of most of the postulated and 
observed behavioral (e.g. adherence) and clinical (e.g. 
depressive symptoms) confounders that may influence 
mortality. Deaths within WIHS are actively ascertained. 
Whereas most prior studies investigating the association 
of race with AIDS death in HAART users have not shown 
a positive association [8—16], in the case of our previous 
WIHS study, we believe we did not meet statistical 
significance as a result of underpowering [21]. 

Study limitations include the potential for unmeasured 
confounders that may also influence AIDS-specific 
clinical outcomes. An additional limitation is the 
inclusion of some self-reported covariates, for example, 
adherence and prior ADI, as well as self-categorized race. 
Self-reported adherence measured on a biannual basis is 
limited by recall bias and the potential for misreporting. 
Additional measures of adherence could include hair 
levels of antiretrovirals [68,69], frequency of filling 
prescriptions, and pill counts, which have their own 
limitations. A previous WIHS study compared self-report 
of AIDS-specific diagnoses with AIDS diagnoses 
documented by county AIDS surveillance registries 
and found a fair degree of accuracy with some variation 
by specific condition [70] . Self-categorization of race may 
not accurately reflect genetic background. Categorization 
of race with ancestry informative markers in WIHS may 
provide a more accurate assessment of race from a genetic 
perspective and is planned for the future. There is also the 
potential for misclassification of death [71—74]. 

In summary, our results provide evidence that African 
American compared to white American women experi- 
ence higher rates of the most important HIV-related 
outcomes, specifically AIDS death and incident ADI. This 
association was specific to AIDS-related outcomes and was 
not explained by many measured potential confounding 
factors including HAART adherence. There is growing 
evidence supporting a differential response to HAART 
possibly based on host genetic variation, including genes 
that influence pharmacokinetics. Future studies replicating 
these associations, and examining specific pharmaco- 
genetic and immuno-genetic variations, as well as 
behaviours, are important and may inform both our 
understanding of HIV pathogenesis and the selection of 
HAART regimens for people of African descent. 
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